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£j \ ABSTRACT 

■ We present a detailed analysis of the two-point correlation function, £(cr, 7r), from the 
' 2dF Galaxy Redshift Survey (2dFGRS). The large size of the catalogue, which contains 

J> , ~ 220 000 redshifts, allows us to make high precision measurements of various properties 

■ of the galaxy clustering pattern. The effective redshift at which our estimates are made is 
rN \ z s w 0.15, and similarly the effective luminosity, L s ps 1.4L*. We estimate the redshift- 

• space correlation function, £(s), from which we measure the redshift-space clustering length, 

so = 6.82 ± 0.28 /i -1 Mpc. We also estimate the projected correlation function, S(cr), and 
the real-space correlation function, £(r), which can be fit by a power-law (r/ro) 1r , with 
r = 5.05 ± 0.26 /i^Mpc, -y r = 1.67 ± 0.03. For r > 20 ft. _1 Mpc, £ drops below a power- 
law as, for instance, is expected in the popular ACDM model. The ratio of amplitudes of the 
real and redshift-space correlation functions on scales of 8 — 30 /i _1 Mpc gives an estimate 
of the redshift-space distortion parameter j3. The quadrupole moment of £ (cr, 7r) on scales 
30 — 40 /i _1 Mpc provides another estimate of /3. We also estimate the distribution function 
of pairwise peculiar velocities, f(v), including rigorously the significant effect due to the in- 
fall velocities, and find that the distribution is well fit by an exponential form. The accuracy 
of our £((7, 7r) measurement is sufficient to constrain a model, which simultaneously fits the 
shape and amplitude of £ (r) and the two redshift-space distortion effects parameterized by (3 
and velocity dispersion, a. We find (3 = 0.49 ± 0.09 and a = 506 ± 52 km s , though the 
best fit values are strongly correlated. We measure the variation of the peculiar velocity dis- 
persion with projected separation, a(a), and find that the shape is consistent with models and 
simulations. This is the first time that (3 and f(v) have been estimated from a self-consistent 
model of galaxy velocities. Using the constraints on bias from recent estimates, and taking 
account of redshift evolution, we conclude that (3(L = L* , z = 0) = 0.47 ± 0.08, and that 
the present day matter density of the Universe, f2 m f=a 0.3, consistent with other 2dFGRS 
estimates and independent analyses. 

Key words: galaxies: statistics, distances and redshifts - large scale structure of Universe - 
cosmological parameters - surveys 
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1 INTRODUCTION 

The galaxy two-point correlation function, f, is a fundamental 
statistic of the galaxy distribution, and is relatively straightforward 
to calculate from observational data. Since the clustering of galax- 
ies is determined by the initial mass fluctuations and their evolu- 
tion, measurements of £ set constraints on the initial mass fluctua- 
tions and their evolution. The astrophysics of galaxy formation in- 
troduces uncertainties, but there is now good evidence that galaxies 
do trace the underlying mass distribution on large scales. 

In this paper we analyse the distribution of ~ 220 000 galaxies 
in the 2-degree Field Galaxy Redshift Survey (2dFGRS, Colless et 
al. 2001). A brief summary of the data is presented in Section 2. 
Much of our error analysis makes use of mock galaxy catalogues 
generated from TV-body simulations which are also discussed in 
Section 2. 

The two-dimensional measurement £(er, 7r), where a is the 
pair separation perpendicular to the line-of-sight and n is the pair 
separation parallel to the line-of-sight provides information about 
the real-space correlation function, the small-scale velocity distri- 
bution, and the systematic gravitational infall into overdense re- 
gions. The spherical average of £(<r, ir) gives an estimate of the 
redshift-space correlation function, £(s), where s = y (n 2 + a 2 ), 
so the galaxy separations are calculated assuming that redshift gives 
a direct measure of distance, ignoring the effects of peculiar ve- 
locities. Integrating £(cr, n) along the line-of-sight sums over any 
peculiar velocity distributions, and so is unaffected by any redshift- 
space effects. The resulting projected correlation function, H(cr), 
is directly related to the real-space correlation function. Our esti- 
mates of 7r), £(s), H(ct) and £(r) are presented in Section 3. 
These statistics have been measured from many smaller redshift 
surveys (e.g. Davis & Peebles 1983; Loveday et al. 1992; Jing, Mo 
& Borner 1998; Hawkins et al. 2001; Zehavi et al. 2002), but since 
they sample smaller volumes, there is a large cosmic variance on 
the results. The large volume sampled by the 2dFGRS leads to sig- 
nificantly more reliable estimates. A preliminary analysis was per- 
formed on the 2dFGRS by Peacock et al. (2001) but we now have 
a far more uniform sample and twice as many galaxies. Madgwick 
et al. (2003) have measured these statistics for spectral-type sub- 
samples of the 2dFGRS. 

Peculiar velocities of galaxies lead to systematic differences 
between redshift-space and real-space measurements, and we can 
consider the effects in terms of a combination of large-scale co- 
herent flows induced by the gravity of large-scale structures, and a 
small-scale random peculiar velocity of each galaxy (e.g. Marzke 
et al. 1995; Jing et al. 1998). The large-scale flows compress the 
contours of £(cr, n) along the n direction, as described by Kaiser 
(1987) and Hamilton (1992). The amplitude of the distortion de- 
pends on the mean density of the universe, fi m , and on how the 
mass distribution is clustered relative to galaxies, which can be pa- 
rameterized in terms of a linear bias b, defined so that S g = b5 m , 
where 8 represents fluctuations in the density field. The random 
component of peculiar velocity for each galaxy means that the ob- 
served £(cr, 7r) is convolved in the 7r coordinate with the pairwise 
distribution of random velocities. In section 4 we describe the con- 
struction of a model £(er, n) from these assumptions about redshift- 
space distortions and also the shape of the correlation function. 

In section 5 we use the Q statistic (Hamilton 1992) based 
on the quadrupole moment of £(<r, tt) to estimate the parameter 
P w Slni 6 /b. In the absence of the small-scale random velocities 
the shape of £(<j, 7r) contours on large scales is directly related to 



the parameter f3. A similar estimate of {3 is provided by the ratio of 
amplitudes of £(s) to f (r) and this is also presented in Section 5. 

In section 6, we use the Landy, Szalay & Broadhurst (1998) 
method to estimate the distribution of peculiar velocities. This tech- 
nique ignores the effect of large-scale distortions and uses the 
Fourier transform of £(<r, 7r) to estimate the distribution of pe- 
culiar velocities, f(v). The large sample volume of the 2dFGRS 
makes our measurements more reliable than previous estimates in 
the same way as for the correlation functions mentioned earlier. 

These two approaches provide reasonable estimates of /3 and 
f(v) so long as the distortions at small and large scales are com- 
pletely decoupled. This is not the case for real data, and so we have 
fitted models which simultaneously include the effects of both /3 
and f(v). The resulting best-fit parameters are the most self con- 
sistent estimates. Previous data-sets have lacked the signal-to-noise 
to allow a reliable multi-parameter fit in this way. Our fitting pro- 
cedure and results are described in Section 7. 

In Section 8, we examine the luminosity and redshift depen- 
dence of /3 and combine our results with estimates of b (Verde et 
al. 2002; Lahav et al. 2002) to estimate S7 m , and compare this with 
other recent analyses. 

In Section 9, we summarise our main conclusions. When 
converting from redshift to distance we assume the Universe 
has a flat geometry with = 0.7, Q m — 0.3 and Hq = 
100 h km s _1 Mpc _1 , so that all scales are in units of fe _1 Mpc. 



2 THE DATA 

2.1 The 2dFGRS data 

The 2dFGRS is selected in the photometric bj band from the APM 
galaxy survey (Maddox, Efstathiou & Sutherland 1990) and its 
subsequent extensions (Maddox et al., in preparation). The bulk 
of the solid angle of the survey is made up of two broad strips, 
one in the South Galactic Pole region (SGP) covering approxi- 
mately -37°.5 < S < -22°5, 21 h 40 m < a < 3 h 40 m and the 
other in the direction of the North Galactic Pole (NGP), spanning 
-7°5 < S < 2°.5, 9 h 50 m < a < 14 h 50 m . In addition to these 
contiguous regions, there are a number of circular 2-degree fields 
scattered randomly over the full extent of the low extinction regions 
of the southern APM galaxy survey. 

The magnitude limit at the start of the survey was set at 
bj = 19.45 but both the photometry of the input catalogue and 
the dust extinction map have been revised since and so there are 
small variations in magnitude limit as a function of position over 
the sky. The effective median magnitude limit, over the area of the 
survey, is bj w 19.3 (Colless et al. 2001). 

The completeness of the survey data varies according to the 
position on the sky because of unobserved fields (mostly around 
the survey edges), un-fibred objects in observed fields (due to col- 
lision constraints or broken fibres) and observed objects with poor 
spectra. The variation in completeness is mapped out using a com- 
pleteness mask (Colless et al. 2001; Norberg et al. 2002a) which is 
shown in Fig.Qfor the data used in this paper. 

We use the data obtained prior to May 2002, which is virtu- 
ally the completed survey. This includes 221 283 unique, reliable 
galaxy redshifts (quality flag ^ 3, Colless et al, 2001). We analyse 
a magnitude-limited sample with redshift limits z m i n = 0.01 and 
Zmax = 0.20, and no redshifts are used from a field with < 70% 
completeness. The median redshift is z mc d ~ 0.11. The random 
fields, which contain nearly 25 000 reliable redshifts are not in- 
cluded in this analysis. After the cuts for redshift, completeness 
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and quality we are left with 165 659 galaxies in total, 95 929 in the 
SGP and 69 730 in the NGP. These data cover an area, weighted 
by the completeness shown in Fig.Q of 647 deg 2 in the SGP and 
446 deg 2 in the NGP, to the magnitude limit of the survey. 

In all of the following analysis we consider the NGP and SGP 
as independent data sets. Treating the NGP and SGP as two inde- 
pendent regions of the sky gives two estimates for each statistic, 
and so provides a good test of the error bars we derive from mock 
catalogues (see below). We have also combined the two measure- 
ments to produce our overall best estimate by simply adding the 
pair counts from the NGP and SGP. The optimal weighting of the 
two estimates depends on the relative volumes surveyed in the NGP 
and SGP, but since these are comparable, a simple sum is close to 
the optimal combination. 

It is important to estimate the effective redshift at which all our 
statistics are calculated. As £ is based on counting pairs of galaxies 
the effective redshift is not the median, but a pair-weighted mea- 
sure. The tail of high redshift galaxies pushes this effective redshift 
to z s w 0.15. Similarly the effective magnitude of the sample we 
analyse is M s — 51og/i ~ —20.0, corresponding to L s w 1.4L* 
(using M* - 5 log h = -19.66, Norberg et al. 2002a). 



2.2 Mock catalogues 

For each of the NGP and SGP regions, 22 mock catalogues were 
generated from the ACDM Hubble Volume simulation (Evrard et 
al. 2002) using the techniques described in Cole et al. (1998), and 
are designed to have a similar clustering signal to the 2dFGRS. 
A summary of the construction methods is presented here but for 
more details see Norberg et al. (2002a) and Baugh et al. (in prepa- 
ration). 

These simulations used an initial dark matter power spectrum 
appropriate to a flat ACDM model with fi m = 0.3 and S1a = 0.7. 
The dark-matter evolution was followed up to the present day and 
then a bias scheme (Model 2 of Cole et al. 1998, with a smooth- 
ing length, Rs = 2 h~ Mpc) was used to identify galaxies from 
the dark matter haloes. The bias scheme used has two free param- 
eters which are adjusted to match the mean slope and amplitude of 
the correlation function on scales greater than a few megaparsec. 
On scales smaller than the smoothing length there is little control 
over the form of the clustering, but in reality the methods employed 
work reasonably well (see later Sections). 

The resulting catalogues have a bias scheme which asymp- 
totes to a constant on large scales, giving f3 = 0.47, but is scale 
dependent on small scales. Apparent magnitudes were assigned to 
the galaxies consistent with their redshift, the assumed Schechter 
luminosity function and the magnitude limit of the survey. The 
Schechter function has essentially the same parameters as in the 
real data (see Norberg et al. 2002a). Then the completeness mask 
and variable apparent magnitude limits were applied to the mock 
catalogues to reproduce catalogues similar to the real data. 

In the analysis which follows, we make use of the real- 
and redshift-space correlation functions from the full Hubble Vol- 
ume simulation. These correlation functions are determined from a 
Fourier transform of the power spectrum of the full Hubble Volume 
cube using the real- and redshift-space positions of the mass parti- 
cles respectively, along with the bias scheme outlined above. This 
allows us to compare our mock catalogue results with that of the 
simulation from which they are drawn to ensure we can reproduce 
the correct parameters. It also allows us to compare and contrast the 
results from the real Universe with a large numerical simulation. 



2.3 Error estimates 

We analyse each of the mock catalogues in the same way as the 
real data, so that we have 22 mock measurements for every mea- 
surement that we make on the real data. The standard deviation 
between the 22 mock measurements gives a robust estimate of the 
uncertainty on the real data. We use this approach to estimate the 
uncertainties for direct measurements from the data, such as the 
individual points in the correlation function, and for best-fit param- 
eters such as so- 

When fitting parameters we use this standard deviation as a 
weight for each data-point and perform a minimum \ 2 analysis 
to obtain the best-fit parameter. The errors that we quote for any 
particular parameter are the rms spread between the 22 best-fit pa- 
rameters obtained in the same way from the mock catalogues. This 
simple way of estimating the uncertainties avoids the complications 
of dealing directly with correlated errors in measured data points, 
while still providing an unbiased estimate of the real uncertainties 
in the data, including the effects of correlated errors. 

Although this approach gives reliable estimates of the uncer- 
tainties, the simple weighting scheme is not necessarily optimal in 
the presence of correlated errors. Nevertheless, for all statistics that 
we consider, we find that the means of the mock estimates agree 
well with the values input to the parent simulations. Also, we have 
applied the technique described by Madgwick et al. (2003) to de- 
correlate the errors for the projected correlation function, using the 
covariance matrix estimated from the mock catalogues. We found 
a 0.1a difference between the best-fit values using the two meth- 
ods. So, we are confident that our measurements and uncertainty 
estimates are robust and unbiased. 



3 ESTIMATES OF THE CORRELATION FUNCTION 

The two point correlation function, £, is measured by comparing 
the actual galaxy distribution to a catalogue of randomly distributed 
galaxies. These randomly distributed galaxies are subject to the 
same redshift, magnitude and mask constraints as the real data and 
we modulate the surface density of points in the random catalogue 
to follow the completeness variations. We count the pairs in bins of 
separation along the line-of-sight, n, and across the line-of-sight, 
a, to estimate £(cr, n). Spherically averaging these pair counts pro- 
vides the redshift-space correlation function £(s). Finally, we es- 
timate the projected function H(cr) by integrating over all velocity 
separations along the line-of-sight and invert it to obtain 

3.1 Constructing a random catalogue 

To reduce shot noise we compare the data with a random catalogue 
containing ten times as many points as the real catalogue. This ran- 
dom catalogue needs to have a smooth selection function matching 
the N(z) of the real data. We use the 2dFGRS luminosity func- 
tion (Norberg et al. 2002a) with M b * T - 5 log h = -19.66 and 
a = —1.21 to generate the selection function, following the change 
in the survey magnitude limit across the sky. When analysing the 
mock catalogues, we use the input luminosity function to generate 
the selection function, and hence random catalogues. 

As an alternative method, we also fitted an analytic form for 
the selection function (Baugh & Efstathiou 1993) to the data, and 
generated random catalogues using that selection function. We have 
calculated all of our statistics using both approaches, and found that 
they gave essentially identical results for the data. When analysing 
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Figure 1. The redshift completeness masks for the NGP (top) and SGP (bottom). The greyscale shows the completeness fraction. 



the mock catalogues, we found that the luminosity function method 
was more robust to the presence of large-scale features in the N(z) 
data. Thus, all of our quoted results are based on random catalogues 
generated using the luminosity function. 

3.2 Fibre collisions 

The design of the 2dF instrument means that fibres cannot be placed 
closer than approximately 30 arcsec (Lewis et al. 2002), and so both 
members of a close pair of galaxies cannot be targeted in a single 
fibre configuration. Fortunately, the arrangement of 2dFGRS tiles 
means that not all close pairs are lost from the survey. Neighbouring 
tiles have significant areas of overlap, and so much of the sky is 
targeted more than once. This allows us to target both galaxies in 
some close pairs. Nevertheless, the survey misses a large fraction of 
close pairs. It is important to assess the impact of this omission on 
the measurement of galaxy clustering and to investigate schemes 
that can compensate for the loss of close pairs. 

To quantify the effect of these so-called 'fibre collisions' we 
have calculated the angular correlation function for galaxies in the 
2dFGRS parent catalogue, w p (0), and for galaxies with redshifts 
used in our £ analysis, w z (8). We used the same mask to determine 
the angular selection and apparent magnitude limit for each sample 
as in Fig.Q Note that the mask is used only to define the area of 
analysis, and the actual redshift completeness values are not used 
in the calculation of w. In our £ analyses we impose redshift lim- 
its 0.01 < z < 0.2, which means that the mean redshift of the 
redshift sample is lower than the parent sample. We used Limber's 
equation (Limber 1954) to calculate the scale factors in amplitude 



<t [h 1 Mpc] 
0.1 1 




e [■] 

Figure 2. Top panel: w(8) for the mean of the NGP and SGP redshift cata- 
logues (solid points), the mean of the masked parent catalogues (solid line), 
and the full APM result (error bars). Bottom panel: The parent catalogue 
result divided by the redshift catalogue results (uncorrected - solid points; 
collision corrected (see Section l33l - open points). The solid line is the 
curve used to correct the fibre collisions. The top axis converts 9 into a 
projected separation, a, at the effective redshift of the survey, z s = 0.15. 
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Figure 3. Redshift distributions, N(z), for the 2dFGRS data (solid lines) and the normalised random catalogues generated using the survey luminosity function 
(dashed lines) for the (a) SGP and (b) NGP. 



and angular scale needed to account for the different redshift dis- 
tributions. The solid line in Fig.|2|shows w p , and the filled circles 
show w z after applying the Limber scale factors. The error bars in 
Fig.|2|show w(6) from the full APM survey (Maddox, Efstathiou 
& Sutherland 1996), also scaled to the magnitude limit of the 2dF- 
GRS parent sample. On scales 8 > 0.03° all three measurements 
are consistent. On smaller scales w z is clearly much lower than 
w p , showing that the fibre collision effect becomes significant and 
cannot be neglected. 

The ratio of galaxy pairs counted in the parent and redshift 
samples is given by (1 + w p )/(l + w z ), which is shown by the 
filled circles in the lower panel of Fig.|2| As discussed in the next 
section, we use this ratio to correct the pair counts in the £ analysis. 



3.3 Weighting 

Each galaxy and random galaxy is given a weighting factor depend- 
ing on its redshift and position on the sky. The redshift dependent 
part of the weight is designed to minimize the variance on the es- 
timated £ (Efstathiou 1988; Loveday et al. 1995), and is given by 
1/(1 + 4irn(zi) J3(s)), where n(z) is the density distribution and 
Js(s) = Jq £(s')s' 2 ds' . We use n(z) from the random catalogue 
to ensure that the weights vary smoothly with redshift. We find that 
our results are insensitive to the precise form of J3 but we derived 
it using a power law £ with so = 13.0 and 7,, = 0.75 and a maxi- 
mum value of J3 = 400. This corresponds to the best-fit power law 
over the range 0.1 < s < 3 h~ 1 Mpc with a cutoff at larger scales. 

We also use the weighting scheme to correct for the galax- 
ies that are not observed due to the fibre collisions. Each galaxy- 
galaxy pair is weighted by the ratio wf — (1 + w p )/(l + w z ) 
at the relevant angular separation according to the curve plotted in 
the bottom panel of Fig.|2| This corrects the observed pair count to 
what would have been counted in the parent catalogue. The open 
points in Fig.|2| which have the collision correction applied, show 
that this method can correctly recover the parent catalogue result 
and hence overcome the fibre collision problem. Since the random 
catalogues do not have any close-pair constraints, only the galaxy- 
galaxy pair count needs correcting in this way. We also tried an 
alternative approach to the fibre-collision correction that we used 
previously in Norberg et al. (2001, 2002b) where the weight for 
each unobserved galaxy was assigned equally to its ten nearest 
neighbours. This produced similar results for 6 > 0.03°, but did 



not help on smaller scales. All of our results are presented using 
the Wf weighting scheme. Hence each galaxy, i, is weighted by the 
factor, 



(1) 



l+47Tn(zi)J 3 (s)' 

and each galaxy-galaxy pair i,j is given a weight WfWiWj, whereas 
each galaxy-random and random-random pair is given a weight 

WiWj. 

3.4 The two-point correlation function, £(cr, n) 
We use the £ estimator of Landy & Szalay (1993), 
DD - 2DR + RR 



RR 



(2) 



where DD is the normalised sum of weights of galaxy-galaxy 
pairs with particular (a, n) separation, RR the normalised sum 
of weights of random-random pairs with the same separation in 
the random catalogue and DR the normalised sum of weights of 
galaxy-random pairs with the same separation. To normalise the 
pair counts we ensure that the sum of weights of the random cat- 
alogue equal the sum of weights of the real galaxy catalogue, as 
a function of scale. We find that other estimators (e.g. Hamilton 
1993) give similar results. 

The N(z) distributions for the data and random catalogues 
(scaled so that the area under the curve is the same as for the ob- 
served data) are shown in Fig. [3] It is clear that N(z) for the ran- 
dom catalogues are a reasonably smooth fit to N(z) for the data. 
Norberg et al. (2002a) showed that large 'spikes' in the N(z) are 
common in the mock catalogues, and so similar features in the data 
redshift distributions indicate normal structure. 

The resulting estimates of £(<j, n) calculated separately for the 
SGP and NGP catalogues are shown in Fig.|4] along with the com- 
bined result. The velocity distortions are clear at both small and 
large scales, and the signal-to-noise ratio is in general very high for 
a and tt values less than 20 h~ 1 Mpc; it is ~ 6 in each 1 h~ 1 Mpc 
bin at s = 20 /i _1 Mpc. At very large separations £(<r, ir) becomes 
very close to zero, showing no evidence for features that could be 
attributed to systematic photometric errors. 

We used an earlier version of the 2dFGRS catalogue to carry 
out a less detailed analysis of £(cr, n) (Peacock et al. 2001). The 
current redshift sample has about 1.4 times as many galaxies, 
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Figure 4. Grey-scale plots of the 2dFGRS §(cr, w) (in 1 h 1 Mpc bins) for (a) the SGP region, (b) the NGP region and (c) the combined data. Contours are 
overlaid at § = 4.0, 2.0, 1.0, 0.5, 0.2 and 0.1. 



though more importantly it is more contiguous, and the revised 
photometry has improved the uniformity of the sample. Neverthe- 
less our new results are very similar to our earlier analysis, demon- 
strating the robustness of our results. The current larger sample al- 
lows us to trace £ out to larger scales with smaller uncertainties. 
Also, in our present analysis we analyse mock catalogues to ob- 
tain error estimates which are more precise than the previous error 
approximation (see Section l7"3l . 



3.5 The redshift-space correlation function, £(s) 

Averaging £(<r, n) at constant s gives the redshift-space correlation 
function, and our results for the NGP and SGP are plotted in Fig. [5] 
on both log and linear scales. The NGP and SGP measurements dif- 
fer by about 2a between 20 and 50 hT 1 Mpc, and we find one mock 
whose NGP and SGP measurements disagree by this much, and so 
it is probably not significant. We tried shifting M* by 0.1 mag to 
better fit the N(z) at z > 0.15 in the SGP, and this moved the data 
points by ~ 0.2a for 20 < s < 50 /i -1 Mpc. 

The redshift-space correlation function for the combined data 
is plotted in Fig.[6jin the top panel. It is clear that the measured £ (s) 
is not at all well represented by a universal power law on all scales, 
but we do make an estimate of the true value of the redshift-space 
correlation length, so, by fitting a localised power-law of the form, 



(3) 



using a least-squares fit to log(£) as a function of log(s), using two 
points either side of f (s) = 1. This also gives a value for the local 
redshift-space slope, 73. The best-fit parameters for the separate 
poles and combined estimates are listed in Table Q In the inset 
of Fig. [6] we can see, at a low amplitude, that £(s) goes negative 
between 50 < s < 90 ft" 1 Mpc. 

In the bottom panel of Fig.|6|we examine the shape of £(s) 
more carefully. The points are the data divided by a small scale 
power law fitted on scales 0.1 < s < 3 ft~ 1 Mpc (dashed line). 
The data are remarkably close to the power-law fit for this lim- 
ited range of scales, and follow a smooth break towards zero for 
3 < s < 60 ft -1 Mpc. The measurements from the Hubble Vol- 



ume simulation are shown by the solid line, and it matches the data 
extremely well on scales s > 4 /i~ 1 Mpc. On smaller scales, where 
the algorithm for placing galaxies in the simulation has little control 
over the clustering amplitude (as discussed in Section l2~2t . there are 
discrepancies of order 50%. 

The mean £(s) determined from the mock catalogues agrees 
well with the true redshift-space correlation function from the full 
Hubble Volume. This provides a good check that our weighting 
scheme and random catalogues have not introduced any biases in 
the analysis. 



3.6 Redshift-space comparisons 

Redshift-space correlation functions have been measured from 
many redshift surveys, but direct comparisons between different 
surveys are not straightforward because galaxy clustering depends 
on the spectral type and luminosity of galaxies (e.g. Guzzo et 
al. 2000; Norberg et al. 2002b; Madgwick et al. 2003). Direct com- 
parisons can be made only between surveys that are based on simi- 
lar galaxy selection criteria. The 2dFGRS is selected using pseudo- 
total magnitudes in the 6j band, and the three most similar surveys 
are the Stromlo-APM survey (SAPM, Loveday et al. 1992), the 
Durham UKST survey (Ratcliffe et al. 1998) and the ESO Slice 
Project (ESP, Guzzo et al. 2000). The Las Campanas Redshift Sur- 
vey (LCRS, Lin et al. 1996, Jing et al. 1998) and Sloan Digital Sky 
Survey (SDSS, Zehavi et al. 2002) are selected in the R band, but 
have a very large number of galaxies, and so are also interesting for 
comparisons. 

The non-power-law shape of £(s) makes it difficult to com- 
pare different measurements of so and 7 S , because the values de- 
pend sensitively on the range of s used in the fitting procedure. In 
Fig. 0a) we compare the £(s) measurements directly for the 2dF- 
GRS, SAPM, Durham UKST and ESP surveys. Our estimate of 
£(s) is close to the mean of previous measurements, but the uncer- 
tainties are much smaller. Although we quote uncertainties that are 
similar in size to previous measurements, we have used the scatter 
between mock catalogues to estimate them, rather than the Pois- 
son or boot-strap estimates that have been used before and which 
seriously underestimate the true uncertainties. 



© 2003 RAS, MNRAS OOO.rnn^l 



The 2dFGRS: correlation functions, peculiar velocities and the matter density of the Universe 7 



* i 



-1 — I — r 



| 20 | | 40 | 50 80 | | 1 00 | | 

<-0A 1 10 100 

s [h MpcJ 

Figure 5. The redshift-space correlation function for the NGP (open points) 
and SGP (solid points) 2dFGRS data with error bars from the rms of mock 
catalogue results. Inset plotted on a linear scale. 
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Figure 6. Top panel: The redshift-space correlation function for the com- 
bined data (points) with error bars from the rms of the mock catalogue re- 
sults. The dashed line is a small scale power law fit, (so = 13 h M.pc, 
7s = 0.75) and the dot-dashed line is the best-fit to points around so. 
(s = 6.82 /i-iMpc, 7 S = 1.57). Inset is on a linear scale. Bottom 
panel: As above, divided by the small scale power law. The solid line shows 
the result from the Hubble Volume simulation. 



Fig. 0b) shows the 2dFGRS measurements together with the 
LCRS and SDSS measurements. On scales s > 4 /i _1 Mpc there 
appears to be no significant differences between the surveys, but 
for s < 2 /i^Mpc the LCRS and SDSS have a higher ampli- 
tude than the 2dFGRS. This difference is likely to be caused by 
the different galaxy selection for the surveys, though the SDSS re- 
sults shown are for the Early Data Release (EDR) and have larger 
errors than the 2dFGRS points. The 2dFGRS is selected using fej, 
whereas the SDSS and LCRS are selected in red bands. Since the 
red (early type) galaxies are more strongly clustered than blue (late 
type) galaxies (e.g. Zehavi et al. 2002; and via spectral type, Nor- 
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Figure 7. Comparison of 2dFGRS f(s) with (a) other fej band selected 
surveys as indicated and (b) R band selected surveys as indicated. These 
results are discussed in the text. 



berg et al. 2002b), we should expect that £ will be higher for red 
selected surveys than a blue selected survey. This issue is examined 
further in Madgwick et al. (2003). 

3.7 The projected correlation function, S(<r) 

The redshift-space correlation function differs significantly from 
the real-space correlation function because of redshift-space dis- 
tortions (see Section 4). We can estimate the real-space correlation 
length, ro, by first calculating the projected correlation function, 
E(<t). This is related to f [a, it) via the equation, 



(a) =2/ t{a,*)d* 



(4) 



though in practice we set the upper limit in this integral to 7r max = 
70 /i _1 Mpc. The result is insensitive to this choice for 7r max > 
60 /i _1 Mpc for our data. Since redshift space distortions move 
galaxy pairs only in the n direction, and the integral represents a 
sum of pairs over all 7r values, S(<r) is independent of redshift- 
space distortions. It is simple to show that E(cr) is directly related 
to the real-space correlation function (Davis & Peebles 1983), 

E(cr) _ 2 [°° r£(r)dr 



(r 2 



(5) 



If the real-space correlation function is a power law this can be 
integrated analytically. We write £(r) = (r/ro)~ Jr , where the 
P superscripts refer to the 'Projected' values, rather than the 'In- 
verted' values which are calculated in Section l3~8l and denoted by 
I. With this notation we obtain, 
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Figure 8. The projected correlation functions for the NGP (open points) 
and SGP (solid points) 2dFGRS data with error bars from the rms spread 
between mock catalogue results. Inset plotted on a linear scale. 




Figure 9. Top panel: The projected correlation function of the combined 
data with error bars from the rms spread between mock catalogue results. 
The dashed line is the best-fit power-law for 0.1 < a < 12 ft" 1 Mpc 
(r = 4.98, 7r = 1.72, A = 3.97). Inset is plotted on a linear scale. 
Bottom panel: The combined data divided by the power-law fit. 



a 



r(|)r( 



7 r r -l- 



r(^) 



(6) 



The parameters 7^ and rjf can then be estimated from the mea- 
sured 3(<r), giving an estimate of the real-space clustering inde- 
pendent of any peculiar motions. 

The projected correlation functions for the NGP and SGP are 
shown in Fig.[8|and the combined data result is shown in Fig.|5| The 
best-fit values of 7^ and ro for 0.1 < a < 12 h~ 1 Mpc are shown 
in Table Q Over this range 3(c) /a is an accurate power law, but 
it steepens for a > 12 h~ 1 Mpc. This deviation from power-law 
behaviour limits the scales that can be probed using this approach. 



Table 1. Best-fit parameters to £. For so an d 7s the fit to £(s) uses only 
points around s = so. For i~q , 7,? and Af/yf 3 ) the fit to 3(<x)/cr uses all 
points with 0.1 < a < 12 h~ 1 Mpc. For and 7^ the fit to the inverted 
£(r) uses all points with 0.1 < r < 12 h~ 1 WLpc. In each case the errors 
quoted are the rms spread in the results obtained from the same analysis 
with the mock catalogues. 



Parameter 


SGP 


NGP 


Combined 


s (/i _1 Mpc) 

7s 


6.92 ± 0.36 
1.51 ± 0.08 


6.72 ±0.41 
1.64 ±0.08 


6.82 ±0.28 
1.57 ±0.07 


(h- 1 Mpc) 


5.05 ± 0.32 
1.68 ± 0.06 
4.17 ±0.23 


4.79 ±0.31 
1.77 ±0.07 
3.77 ±0.28 


4.95 ±0.25 
1.72 ±0.04 
3.99 ±0.16 


rl (/i _1 Mpc) 


5.09 ±0.35 
1.65 ± 0.03 


5.08 ±0.28 
1.70 ±0.04 


5.05 ± 0.26 
1.67 ±0.03 



3.8 The real-space correlation function, £(r) 

It is possible to estimate f (r) by directly inverting H(cr) without 
making the assumption that it is a power law (Saunders, Rowan- 
Robinson & Lawrence 1992, hereafter S92). They recast Eqn. [5] 
into the form, 



50 {dE{a)/da) 
(a 2 — r 2 ) 2 



da. 



(7) 



Assuming a step function for H(cr) = Hi in bins centered on <j;, 
and interpolating between values, 



7T ^ ' 



In 



j+i 



7r 1 — ' a-j 



(8) 



for r = <7j. S92 suggest that their method is only good for scales 
r < 30 /i- 1 Mpc in the QDOT survey because r becomes compa- 
rable to the maximum scale out to which they can estimate S. We 
can test the reliability of our inversion of the 2dFGRS data using 
the mock catalogues. 

In Fig. E3 we show the mean £(r) as determined from the 
mock catalogues using the method of S92. We compare this to 
the real-space correlation function determined directly from the 
Hubble Volume simulation, from which the mock catalogues are 
drawn. The agreement is excellent and shows that the method 
works and that we can recover the real-space correlation function 
out to 30 h~ Mpc. Like S92 we find that beyond this scale the 
method begins to fail and the true £(r) is not recovered. 

We have applied this technique to the combined 2dFGRS data 
and obtain the real-space correlation function shown in Fig. II II 
The data are plotted out to only 30 /i^Mpc due to the limi- 
tations in the method described above. On small scales £(r) is 
well represented by a power law, and a best-fit over the range 
0.1 < r < 12 /i _1 Mpc gives the results for r 7 ] and 7^ shown 
in TableQ 

The points in the bottom panel of Fig.^3 snow me 2dFGRS 
data divided by the best-fit power law. It can be seen that at scales 
0.1 < r < 20 h~ 1 Mpc the data £(r) is close to the best-fit power- 
law but does show hints of non power-law behaviour (see also dis- 
cussion below). 
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r [h Mpc] 

Figure 10. The mean real-space correlation function determined from the 
22 mock catalogues using the method of S92. The solid line is the true f (r) 
from the Hubble Volume and the agreement is excellent. Note the changed 
scale from previous plots. 




r [h 'Mpc] 

Figure 11. Top panel: The real-space correlation function of the com- 
bined 2dFGRS using the method of S92, with error bars from the rms 
spread between mock catalogues. The dashed line is the best-fit power law 
(r"o = 5.05, 7 r = 1.67). Inset is plotted on a linear scale. Bottom panel: 
The data divided by the power law fit. The solid line is the de-projected 
APM result (Padilla & Baugh 2003) as discussed in the text. The dotted 
line is the result from the Hubble Volume. 



3.9 Real-space comparisons 

In the inverted £(r) (and possibly H[er]) there is a weak excess of 
clustering over the power-law for 5 < r < 20 /i _1 Mpc. This 
has been previously called a 'shoulder' in £ (see e.g. Ratcliffe et 
al. 1998). Though the amplitude of the feature in our data is rather 
low, it has been consistently seen in different surveys, and proba- 
bly is a real feature. After submission of this work, Zehavi et al. 
(2003) also saw this effect in the SDSS projected correlation func- 



Table 2. Measurements of §(r) from 2dFGRS and other surveys, with the 
quoted uncertainties as published. The various authors have used very dif- 
ferent ways to estimate errors though none have included the effects of cos- 
mic variance. Since we have used the scatter between mock catalogues to 
take account of this, our estimates are actually dominated by cosmic vari- 
ance. These results are measured at different effective luminosities, red- 
shifts and for different galaxy types. 



Survey 



r (fc-iMpc) 7? . 



2dFGRS (P) 


4.95 ±0.25 


1.72 ±0.04 


2dFGRS (/) 


5.05 ±0.26 


1.67 ±0.03 


SAPM 


5.1 ±0.3 


1.71 ±0.05 


ESP 


4.15 ±0.2 


1 67+ 007 

1 - D '-0.09 


Durham UKST 


5.1 ±0.3 


1.6 ±0.1 


LCRS 


5.06 ±0.12 


1.86 ±0.03 


SDSS 


6.14 ±0.18 


1.75 ±0.03 



tion and explained the inflection point as the transition scale be- 
tween a regime dominated by galaxy pairs in the same halo and 
a regime dominated by pairs in separate haloes. Magliochetti & 
Porciani (2003) have found the same effect when examining corre- 
lation functions of different types of 2dFGRS galaxy. 

The dotted line in the bottom panel of Fig. 1111 shows the 
Hubble Volume simulation which agrees well with the data for 
?• > 1 h~ 1 Mpc. On smaller scales the Hubble Volume £ shows 
significant deviations from a power-law. On these scales, the galaxy 
clustering amplitude in the simulation is incorrectly modelled since 
the assignment of galaxies to particles is based on the mass distri- 
bution smoothed on a scale of 2 ft _1 Mpc (as discussed in Sec- 
tion |2j2j- The solid line in the bottom panel of Fig. 1 1 1 1 is the de- 
projected APM result (Padilla & Baugh 2003), scaled down by a 
factor (1 + z s ) a , with a — 1.7, suitable for evolution in a ACDM 
cosmology. There is good agreement between the 2dFGRS and 
APM results which are obtained using quite different methods. 

We have estimated ro and y r by fitting to the projected cor- 
relation function E(<t)/<j, and also by inverting S(cr) /a and then 
fitting to £(r). The best-fit values from the two methods are shown 
in Table0 and it is clear they lead to very similar estimates of ro 
and 7 r . This confirms that the power-law assumption in Section l3~71 
is a good approximation over the scales we consider. 

Table|2|lists ro and 7,. for the 2dFGRS and other surveys es- 
timated using power-law fits to the projected correlation function 
H(cr). As mentioned in Section|^|the SAPM, Durham UKST and 
the ESP are bj selected surveys, and so should be directly compa- 
rable to the 2dFGRS. The values of ro and j r for these surveys 
all agree to within one standard deviation, except ro for the ESP, 
which appears to be significantly lower. It is likely that the quoted 
uncertainties for the ESP and Durham UKST parameters are under- 
estimated since they did not include the effect of cosmic variance. 
Since they each sample relatively small volumes, this will be a large 
effect. The sparse sampling strategy used in the SAPM means that 
it has a large effective volume, and so the cosmic variance is small. 

As in Section l3~5l the red-selected surveys, LCRS and SDSS, 
are significantly different from the other surveys. The discrepancies 
are most likely due to the fact that the amplitude of galaxy cluster- 
ing depends on galaxy type, and that red-selected surveys have a 
different mix of galaxy types. We can make a very rough approx- 
imation of the expected change in £ by considering how the mean 
colour difference of early and late populations changes the relative 
fraction of the two populations when a magnitude limited sample 
is selected in different pass bands. Zehavi et al. (2002) split their 
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r— selected SDSS sample into 19603 early-type galaxies and 9532 
late-type galaxies. The mean (g — r) colours are 0.5 and 0.9 respec- 
tively. The 2dFGRS is selected using bj which is close to g, and so 
compared to the r selection, the median depth for blue galaxies will 
be larger that for red galaxies. The number of early and late types 
will roughly scale in proportion to the volumes sampled, and so 
the ratio of early-to-late galaxies in the 2dFGRS will be roughly 
~ (19603/9532) x io ^ 5 " 9 ) = 1.18. Note that this colour 
split leads to a very different ratio of early-to-late galaxies com- 
pared to the split used by Madgwick et al. (2003). Assuming the 
early and late correlation functions trace the same underlying field, 
the combined correlation function will be 



■ Nearly Nearly ~\~ ^late^late i 

^ tot ~ I 7, ~Z7, ' ^ n 

Nearly T dilate 



(9) 



From the power law fits of Zehavi et al., the ratio of bias val- 
ues at 1 Mpc is &eariy/&iate = 4.95. Inserting the different ratios 
nearly /niato appropriate to the red and blue selected samples we 
find that the expected ratio of £ for a red selected sample com- 
pared to a blue selected sample is roughly 1.36. Scaling the 2dF- 
GRS values of r = 5.05 and j r = 1.67 leads to a SDSS value 
of ro = 5.95 for 7 r = 1.75, within la of the actual SDSS value. 
This simple argument indicates that the observed difference in £ 
between the red and blue selected surveys is consistent with the 
different population mixes expected in the surveys. The extra sur- 
face brightness selection applied to the LCRS may also introduce 
significant biases. 

Each survey is also likely to have a different effective lumi- 
nosity and, as has been shown by Norberg et al. (2001), this will 
cause clustering measurements to differ. The relation for 2dFGRS 
galaxies found by Norberg et al. (2001) was, 



ro 



0.85 + 0.15 



(10) 



which gives, for L = 1AL* (see Section^}, rg = 4.71 ± 0.24, 
which will allow direct comparisons with other surveys. 



4 REDSHIFT-SPACE DISTORTIONS 

When analysing redshift surveys it must be remembered that the 
distance to each galaxy is estimated from its redshift and is not the 
true distance. Each galaxy has, superimposed on its Hubble motion, 
a peculiar velocity due to the gravitational potential in its local en- 
vironment. These peculiar velocities can be in any direction and, 
since this effect distorts the correlation function, it can be used to 
measure two important parameters. 

The peculiar velocities are caused by two effects. On small 
scales, random motions of the galaxies within groups cause a ra- 
dial smearing known as the 'Finger of God'. On large scales gravi- 
tational instability leads to coherent infall into overdense regions 
and outflow from underdense regions. We analyse the observed 
redshift-space distortions by modeling £(a,ir). We start with a 
model of the real-space correlation function, £ (r), and include the 
effects of large-scale coherent infall, which is parameterized by 
/3 ~ fi„; 6 /b, where b is the linear bias parameter. We then con- 
volve this with the form of the random pairwise motions. 



4.1 Constructing the model 

Kaiser (1987) pointed out that, in the linear regime, the coherent in- 
fall velocities take a simple form in Fourier space. Hamilton (1992) 
translated these results into real space, 



C'(a,^)=£o(s)Po(/x)+e2( S )P2(M)+e4( S )P4(/x) 



(ID 



where Pf(/i) are Legendre polynomials, = cos(#) and 8 is the 
angle between r and n. The relations between £e, £(r) and (3 for a 
simple power-law £(r) = (r/ro) ,r are (Hamilton 1992), 



1 + f + y ) 



'4/3 4/3 2 



7r 



Jr - 3 

7 r (2 + 7 r ) 
35 V (3 - 7r)(5 - 1r) 



CM 

£(r). 



(12) 



(13) 



(14) 



The Appendix has more details of this derivation and gives the 
equations for the case of non-power law forms of £. 

We use these relations to create a model £'(<x, it) which we 
then convolve with the distribution function of random pairwise 
motions, f(v), to give the final model £(er, n) (Peebles 1980): 



/oo 
£'((j, n — v/H )f(v)dv 
-oo 



(15) 



and we choose to represent the random motions by an exponential 
form, 



tr \ 1 ( 



(16) 



where a is the pairwise peculiar velocity dispersion (often known as 
(T12). An exponential form for the random motions has been found 
to fit the observed data better than other functional forms (e.g. Rat- 
cliffe et al. 1998; Landy 2002; see also Section|6j. 



4.2 Model assumptions 

In this model we make several assumptions. Firstly, we assume a 
power-law for the correlation function. The power-law approxima- 
tion is a good fit on scales < 20 h~ 1 Mpc but is not so good at 
larger scales. This limits the scales which we can probe using this 
method. In SectionQ we consider non-power-law models for £(r), 
and recalculate EQns. [l2l to ll4l using numerical integrals (see Ap- 
pendix), allowing us to reliably use scales > 20 /i _1 Mpc. Sec- 
ondly, we assume that the linear theory model described above 
holds on scales < 8 h~ Mpc, which is almost certainly not true. 
We also consider this in Section 7. Finally, we assume an expo- 
nential distribution of peculiar velocities with a constant velocity 
dispersion, a, (Ean. ll6> and this is discussed and justified in Sec- 
tion[6]and Sectionl7~4l 



4.3 Model plots 

To illustrate the effect of redshift-space distortions on the £(cr, 7r) 
plot we show four model £(<r, 7r)'s in Fig. 1121 If there were no dis- 
tortions, then the contours shown would be circular, as in the top 
left panel due to the isotropy of the real-space correlation function. 
On small a scales the random peculiar velocities cause an elonga- 
tion of the contours in the n direction (the bottom left panel). On 
larger scales there is the flattening of the contours (top right panel) 
due to the coherent infall. The bottom right panel is a model with 
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Figure 12. Plot of model £ (<r, tt) 's calculated as described in Sectionl4l The 
lines represent contours of constant £(cr, it) = 4.0, 2.0, 1.0, 0.5, 0.2 and 0.1 
for different models. The top left panel represents an undistorted correlation 
function (a = 0, f) = 0), the top right panel is a model with coherent infall 
added (a = 0, (3 = 0.4), the bottom left panel is a model with just random 
pairwise velocities added (a = 500 km s _1 , (3 = 0) and the bottom right 
panel has both infall and random motions added (a = 500 km s^ 1 , (3 = 
0.4). These four models have rg = 5.0 ft - 1 Mpc and -y r = 1.7. 



both distortion effects included. Comparing the models of £(a, n) 
to the 2dFGRS measurements in Fig.[4]it is clear that the data show 
the two distortion effects included in the models. In Section 0wc 
use the data to constrain the model directly, and deduce the best-fit 
model parameters. 



5 ESTIMATING (3 

Before using the model described above to measure the parameters 
simultaneously, we first use methods that have been used in pre- 
vious studies. This allows a direct comparison between our results 
and previous work. 



5.1 Ratio of fs 

The ratio of the redshift-space correlation function, £(s), to the 
real-space correlation function, £(r), in the linear regime gives an 
estimate of the redshift distortion parameter, /3 (see Eqn, ll2t . 

frH + f-4- (17) 

f(r) 3 5 

Our results for the combined 2dFGRS data, using the inverted form 
of £(r), are shown in Fig. ll3l bv the solid points. The mean of the 
mock catalogue results is shown by the white line, with the rms 
errors shaded and the estimate from the Hubble Volume is shown 
by the solid line. The data are consistent with a constant value, and 
hence linear theory, on scales > 4 /i -1 Mpc. 

The mock catalogues and Hubble Volume results asymptote to 
(3 = 0.47, the true value of (3 in the mocks. The 2dFGRS data in 
the range 8 — 30 /i Mpc are best-fit by a ratio of 1.34 ± 0.13, 




Figure 13. The ratio of f (s) to £(r) for the 2dFGRS combined data (solid 
points), and the Hubble Volume (solid line). The mean of the mock cata- 
logue results is also shown (white line), with the rms errors shaded. The 
error bars on the 2dFGRS data are from the rms spread in mock catalogue 
results. 



corresponding to (3 — 0.45 ± 0.14. The maximum scale that we 
can use in this analysis is determined by the uncertainty on £(r) 
from the Saunders et al. inversion method discussed in Section lX8l 



5.2 The quadrupole moment of £ 

We now measure (3 using the quadrupole moment of the correlation 
function (Hamilton 1992), 



Q(s) = 



1/3+ f/3 2 
l + f/34- 



X fo(s')s' 2 ds'-£o(* 



where £e is given by, 
, . . 21+1 

&00 = — ;t— 



£,(a,n)Pi(fi)dfi. 



(18) 



(19) 



These equations assume the random peculiar velocities are negli- 
gible and hence measuring Q gives an estimate of (3. The random 
uncertainties in this method are small enough that we obtain reli- 
able estimates on scales < 40 h~ Mpc, as shown by the mock 
catalogues (see below), but the data are noisy beyond these scales. 

Fig. 1141 shows Q estimates for the combined 2dFGRS data 
with the inset showing the NGP and SGP separately. The effect of 
the random peculiar velocities can be clearly seen at small scales, 
causing Q to be negative. The best-fit value to the combined data 
for 30 - 40 /i _1 Mpc is Q = 0.55 ± 0.18, which gives a value for 
P — 0.47lg'ig, where the error is from the rms spread in the mock 
catalogue results. The solid line represents a model with (3 = 0.49 
and a — 506 km s" 1 , which matches the data well (see Section 
17. U . Although asymptoting to a constant, the value of Q in the 
model is still increasing at 40 /i _1 Mpc. This shows that non-linear 
effects do introduce a small systematic error even at these scales 
though this bias is small compared to the random error. 

To check whether this method can correctly determine (3 we 
use the mock catalogues. The data points in Fig. ll5l are the mean 
values of Q from the mock catalogues, with error bars on the mean, 
and the dashed line is the true value of S3 = 0.47. The data points 
seem to converge on large scales to the correct value of Q. Fitting 
to each mock catalogue in turn for 30 — 40 /i _1 Mpc gives a mean 
Q = 0.51 + 0.18, corresponding to [3 



0.4315;^. As the models 



© 2003 RAS, MNRAS 000,[T][T9] 



12 Hawkins et al. (The 2dFGRS Team) 




10 20 30 40 



Scale [h 'Mpc] 

Figure 14. The Q factor for the combined 2dFGRS data, restricted to scales 
not dominated by noise with error bars from the rms spread in mock cata- 
logue results. The two dashed lines show the expected answer for different 
values of /3 which approximate the lcr errors. The solid line shows a model 
with (3 = 0.49 and a = 506 km s — 1 (see Section I7T1 . The inset shows 
the result for the NGP (solid points) and SGP (open points); the error bars 
are placed alternately to avoid confusion. 
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Figure 15. The mean Q factor for the mock catalogues with error bars from 
the rms spread in mock catalogue results. The dashed line is the true value 
of f3 = 0.47. The slight bias is caused by the random peculiar velocities as 
discussed in the text. 

showed, the random velocities will lead to an underestimate of j3 
even at 40 /i _1 Mpc, causing the difference between the measured 
and true values. This all shows that we can determine [3 with a slight 
bias but the error bars are large compared to the bias. 

The Q estimates from the individual mock catalogues show 
a high degree of correlation between points on varying scales and 
so the overall uncertainty in Q from averaging over all scales > 
30 h~ 1 Mpc is not much smaller than the uncertainty from a single 
point. It is this fact which makes the spread in results from the 
mock catalogues vital in the estimation of the errors on our result 
(also see Section fOl . 

6 THE PECULIAR VELOCITY DISTRIBUTION 

To this point we have assumed that the random peculiar velocity 
distribution has an exponential form (Eqn. ll6> . This form has been 



used by many authors in the past and has been found to fit the data 
better than other forms (e.g. Ratcliffe et al. 1998). We test this for 
the 2dFGRS data by following a method similar to that of Landy, 
Szalay and Broadhurst (1998, hereafter LSB98). To extract the pe- 
culiar velocity distribution, we need to deconvolve the real-space 
correlation function from the peculiar velocity distribution. 



6.1 The method 

We first take the 2-d Fourier transform of the £(cr, n) grid to give 
^(kajk-rr) and then take cuts along the k a and axes which we de- 
note by T,(k) and II(fc) respectively, so = £(fc CT = k, k n = 
0) and fl(fc) = £(kcr — 0, k n = k). By the slicing-projection 
theorem (see LSB98) these cuts are equivalent to the Fourier trans- 
forms of the real-space projections of £(<r, tt) onto the a and n 
axes. The projection of £ (a, n) onto the a axis is a distortion free 
measurement of 5 but the projection onto the tt axis gives us S 
convolved with the peculiar velocity distribution, ignoring the ef- 
fects of large-scale bulk flows. Since a convolution in real space 
is a multiplication in Fourier space, the ratio of £(fc) to H(k) is 
the Fourier transform, of the velocity distribution that we 

want to estimate. All that is left is to inverse Fourier transform this 
ratio to obtain the peculiar velocity distribution, f(v). LSB98 cut 
their dataset at 32 /i _1 Mpc and applied a Hann smoothing win- 
dow; we use all the raw data. Landy (2002, hereafter L02) used the 
LSB98 method on the 100k 2dFGRS Public Release data and his 
results are discussed below. 

Fitting an exponential to the resulting f(v) curve gives a value 
for a assuming that the infall contribution to the velocity distribu- 
tion is negligible. LSB98 and L02 claim that their method is not 
sensitive to the infall velocities. We show here that this is not the 
case. The additional structure in the Fourier transform of the veloc- 
ity distribution found by L02 is a direct consequence of the infall 
velocities. 



6.2 Testing the models 

To test the LSB98 method we apply the technique to our models, 
described in Section[4] with and without a /3 = 0.4 infall factor, us- 
ing various scales, and with and without a Hann window. In Fig. 1161 
we show the Fourier transform of the peculiar velocity distribution 
and in Fig. ll7l we show the distribution function itself. 

It is clear from Fig. 1161 that the shape of the Fourier trans- 
form at small k is quite badly distorted by the infall velocities. This 
leads to a systematic error in the actual velocity distribution as seen 
in Fig.^| where the measured peculiar velocity dispersions are bi- 
ased low, especially in the case where a smoothing window and a 
limited range of scales are used. In particular the peak of the Fourier 
transform is not at k — 0, and the inferred f(v) goes negative for 
a range of velocities (dashed lines in the lower panels of Fig. ll7> . 
This clearly cannot be interpreted as a physical velocity distribu- 
tion; the method infers negative values because the input model 
£(<j, 7r) is not consistent with the initial assumption of the method, 
which is that all of the distortion in £(er, 7t) is due to random pe- 
culiar velocities. We conclude that both types of peculiar velocity 
need to be considered when making these measurements, and so 
our preferred results come from directly fitting to £(<r, 7r). 

A further complication with the real data is that f(v) may de- 
pend on the pair separation (see discussion in Section 7.3). The 
solid line in Fig. ll8l shows T\f{v)] for a model where a varies from 
500 km s _1 at a = to 300 km s _1 at a = 20 /i _1 Mpc. This is 
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Figure 16. The Fourier transform of the peculiar velocity distribution for 
various parameters (see labels). The solid line is for a model with no 
smoothing and using all scales < 70 ft Mpe, The dashed line is for a 
model cut at 32 h _1 Mpc and smoothed with a Hann window (like L02). 
The dotted line is the Lorentzian equivalent of the input exponential peculiar 
velocity distribution (this coincides with the solid line in the top panels). 
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Figure 17. The peculiar velocity distribution for various models (see la- 
bels). The solid line is the recovered distribution for a model with no 
smoothing and using all scales < 70 h~ 1 Mpc. The dashed line is the 
recovered distribution for a model cut at 32 h _1 Mpc and smoothed with 
a Hann window. The dotted line is the input exponential peculiar velocity 
distribution (this coincides with the solid line in the top panels). 



compared to a model with a = 500 km s (dashed line), a model 
with a = 300 km s -1 (dotted line) and a model where a varies 
from 300 km s" 1 at a = to 500 km s" 1 at a = 20 /i _1 Mpc 
(dot-dashed line). The models with varying a are very close to 
their respective constant a models at all k values, showing that this 
method leads to an estimate of T[f(y)] determined mainly by the 
value of a at small a. 



6.3 The mock catalogues 

The mean of the peculiar velocity distributions for the mock cata- 
logues is shown in Fie. 1191 The distribution is compared to a model, 
shown as the solid line, with j3 — 0.47 and an exponential f(v), 
with dispersion, a = 575 km s - . The exact form of the peculiar 
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Figure 18. The Fourier transform of the peculiar velocity distribution for 
a model with a = 500 km s —1 (dashed line), a = 300 km s _1 (dot- 
ted line). Also shown is a model with a decreasing from 500 km s -1 
to 300 km s _1 from a = to a = 20 h _1 Mpc (solid line) and a 
model with a increasing from 300 km s~ 1 to 500 km s — 1 from a = to 
cr = 20 fe _1 Mpc (dot-dashed line). /3 = 0.5 for all four models. 
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Figure 19. The recovered velocity distribution for the mock catalogues. 
Filled points are the mean result with error bars from the scatter between 
catalogues. This is compared to a pure exponential distribution with a = 
575 km s _1 (dashed line) and a model with a = 575 km s — 1 and (3 = 
0.47 (solid line). 



velocities in the Hubble Volume, and hence mock catalogues, is not 
explicitly specified and it should not be expected to conform to this 
model exactly. 



6.4 The 2dFGRS data 

The Fourier transform of the peculiar velocity distribution for the 
combined 2dFGRS data are shown in Fig. l20l compared to a best-fit 
model with (3 = 0.49 ± 0.05 and a = 570 ± 25 km s _1 . Fie.lzTl 
shows the peculiar velocity distribution itself compared to the same 
model. We showed in Section l6"2l (wifh Fig. 1181 that this was likely 
to be the value of a at small a. The distribution of random pairwise 
velocities does appear to have an exponential form, with a j3 influ- 
ence. Sheth (1996) and Diaferio & Geller (1996) have shown that 
an exponential peculiar velocity distribution is a result of gravita- 
tional processes. 

Ignoring the infall L02 found a — 331 km s" 1 , using the 
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Figure 20. The Fourier transform of the peculiar velocity distribution for 
the combined 2dFGRS data (solid points) compared to a model distribution 
with a = 570 km s _1 and f3 = 0.49 (solid line). 
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Figure 21. The combined 2dFGRS peculiar velocity distribution (solid 
points), compared to a pure exponential distribution with a = 570 km s _1 
(dashed line) and a model with a = 570 km s _1 and /3 = 0.49 (solid 
line). The error bars are from the scatter in mock results. 



smaller, publicly available, sample of 2dFGRS galaxies. We made 
the same approximations and repeated his procedure on our larger 
sample, and find a = 370 km s - . Using our data grid out to 
70 h~ 1 Mpc, with no smoothing and ignoring /3, gives a = 
457 km s -1 . We have shown that the result in L02 is biased 
low by ignoring (3 and that the infall must be properly consid- 
ered in these analyses. As shown in Fig. 1211 our data are reason- 
ably well described by an exponential model with (3 = 0.49 and 
a = 570 km s -1 . 



7 FITTING TO THE £ (er, tt ) GRID 
7.1 Results 

We now fit our £(<j, n) data grid to the models described in Section 
[4] assuming a power-law form for the real-space correlation func- 
tion. This model has four free parameters, (3, ro, 7r and a. The fits 
to the data are done by minimising 



for s < 20 h~ 1 Mpc, where St; is the rms of £ from the mock 
catalogues for a particular a and 7r. This is like a simple \ 2 min- 
imization, but the points are not independent. We tried a fit to £ 
directly but found that it gave too much weight to the central re- 
gions and so instead we fit to log[l + £] so that the overall shape 
of the contours has an increased influence on the fit. The best-fit 
model parameters are listed in Table|3| The errors we quote are the 
rms spread in errors from fitting each mock catalogue in the same 
way. 

There are two key assumptions made in the construction of 
these models. Firstly, although the contours match well at small 
scales, there are good reasons to believe that our linear theory 
model will not hold in the non-linear regime for s < 8 /i _1 Mpc. 
Secondly, we have assumed the power-law model for £(r) and we 
have seen evidence that this is not completely realistic. Using non- 
power law forms will also allow us to probe to larger scales. 

To test whether our result is robust to these assumptions 
we firstly reject the non-linear regime corresponding to s < 
8 /i _1 Mpc. Then, we use the shape of the Hubble Volume £(r) 
instead of a power-law, and finally we extend the maximum scale 
to s = 30 /i _1 Mpc. We showed in Section ^1 that the Hub- 
ble Volume shape gives a good match to the data over the range 
8 < s < 30 ft _1 Mpc (the Appendix gives the relevant equa- 
tions for performing the (3 infall calculation without a power-law 
assumption). 

We find that the best-fit parameters change very little with 
these changes but when using the Hubble Volume £(r), the quality 
of the fit improves significantly. The best-fit model is compared to 
the data in Fig. [22] Notice the excellent agreement on small scales 
even though they are ignored in the fitting process. The best-fit pa- 
rameters are listed in Table|3| and we adopt these results as our final 
best estimates finding (3 = 0.49 ± 0.09. 

If we repeat our analysis on the mock catalogues we find a 
mean value of j3 = 0.475 ± 0.090 (cf. the expected value of 
f3 = 0.47, Section l2"2l . showing that we can correctly determine (3 
using this type of fit. When fitting the mock catalogues it became 
clear that (3 and a are correlated in this fitting procedure, as we 
have seen already with other methods. We use the mock catalogues 
to measure the linear correlation coefficient, r (Press et al. 1992), 
which quantifies this correlation, and find that, between (3 and a, 
r = 0.66. If we knew either parameter exactly, the error on the 
other would be smaller than quoted. 

We also tried other analytical forms for the correlation func- 
tion and also different scale limits and found that some combina- 
tions shifted the results by ~ ler. 



7.2 Comparison of methods 

We have now estimated the real-space clustering parameters using 
three different methods. In Section lT9l we saw that the projection 
and inversion methods gave essentially identical results for ro and 
7 r whereas using 2-d fits we get slightly higher values for ro. 

If ^ (r) was a perfect power-law the different methods would 
give unbiased results for the parameters, but we have seen evi- 
dence that this assumption is not true. The methods therefore, give 
different answers as a result of the different scales and weight- 
ing schemes used, as well as the vastly different treatments of the 
redshift-space distortions. 
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Figure 22. Contours of £ (<r, it) for the 2dFGRS combined data (solid lines) 
and the best-fit model (see Tablel3l using the Hubble Volume £(r) fitted to 
scales 8 < s < 30 /i" 1 Mpc (dashed lines). Contour levels are at £ = 4.0, 
2.0, 1.0, 0.5, 0.2, 0.1, 0.05 and 0.0 (thick line). 



Table 3. Best-fit parameters to the £(er, 7r) grids with errors from the rms 
spread in mock catalogue results. 



Parameter SGP NOP Combined 

Power law £(r): (0 < s < 20 h _1 Mpc) 

(3 0.53 ±0.06 0.48 ±0.08 0.51 ± 0.05 

r (ft-'Mpc) 5.63 ±0.26 5.52 ±0.29 5.58 ±0.19 

7 r 1.66 ±0.06 1.76 ±0.07 1.72 ±0.05 

a(kms- 1 ) 497 ±24 543 ± 26 522 ± 16 

Power law £(r): (8 < s < 20 h^Mpc) 

(3 0.45 ±0.10 0.35 ±0.12 0.49 ± 0.09 

r (h-'Mpc) 6.03 ±0.36 6.06 ± 0.41 5.80 ± 0.25 

7 r 1.74 ±0.08 1.88 ±0.10 1.78 ±0.06 

a(kms- 1 ) 457 ±49 451 ± 51 514 ± 31 

Hubble Volume £(r): (8 < s < 20 h^Mpc) 

(3 0.47 ±0.12 0.50 ±0.14 0.49 ±0.10 

a(kms- 1 ) 446 ±73 544 ± 67 495 ± 46 

Hubble Volume £{r): (8 < s < 30 h~ 1 Mpc) 

13 0.48 ±0.11 0.47 ±0.13 0.49 ± 0.09 

atkms- 1 ) 450 ± 81 545 ± 85 506 ± 52 



7.3 Previous 2dFGRS results 

It is worth contrasting our present results with those obtained in a 
previous 2dFGRS analysis (Peacock et al. 2001). This was based 
on the data available up to the end of 2000: a total of 141 402 
redshifts. The chosen redshift limit was z max = 0.25, yielding 
127 081 galaxies for the analysis of £(<r, n). The present analysis 
uses 165 659 galaxies, but to a maximum redshift of 0.2. Because 
galaxies are given a redshift-dependent weight, this difference in 



redshift limit has a substantial effect on the volume sampled. For 
a given area of sky, changing the redshift limit from z max = 0.2 
to Zmax = 0.25 changes the total number of galaxies by a factor 
of only 1.08, whereas the total comoving volume within z max in- 
creases by a factor of 2. Allowing for the redshift-dependent weight 
used in practice, the difference in effective comoving volume for a 
given area of sky due to the variation in redshift limits becomes a 
factor of 1.6. Since the effective area covered by the present data 
is greater by a factor of 165 659/(127 081/1.08) = 1.4, the total 
effective comoving volume probed in the current analysis is in fact 
15% smaller than in the 2001 analysis; this would suggest random 
errors on clustering statistics about 7% larger than previously. Of 
course, the lower redshift limit has several important advantages: 
uncertainties in the selection function in the tail of the luminos- 
ity function are not an issue (see Norberg et al. 2002a); also, the 
mean epoch of measurement is closer to z = 0. Given that the sky 
coverage is now more uniform, and that the survey mask and selec- 
tion function have been studied in greater detail, the present results 
should be much more robust. 

The other main difference between the present work and that 
of Peacock et al. (2001) lies in the method of analysis. The ear- 
lier work quantified the flattening of the contours of £(cr, n) via the 
quadrupole-to-monopole ratio, £,2(s)/£o(s). This is not to be con- 
fused with the quantity Q(s) from Section l5~2l which uses an inte- 
grated clustering measure instead of £o(s). This is inevitably more 
noisy, as reflected in the error bar, 5/3 = 0.17, resulting from that 
method. The disadvantage of using £2(s)/£o(s) directly, however, 
is that the ratio depends on the true shape of £(r). In Peacock et al. 
(2001), this was assumed to be known from the de-projection of an- 
gular clustering in the APM survey (Baugh & Efstathiou 1993); in 
the present paper we have made a detailed internal estimate of £(r), 
and considered the effect of uncertainties in this quantity. Apart 
from this difference, the previous method of fitting to £2(s)/£o(s) 
should, in principle, give results that are similar to our full fit to 
£(<j, 7r) in Section lTTI The key issue in both cases is the treatment 
of the errors, which are estimated in a fully realistic fashion in the 
present paper using mock samples. The previous analysis used two 
simpler methods: an empirical error on £2(s)/£o(s) was deduced 
from the NGP-SGP difference, and correlated data were allowed 
for by estimating the true number of degrees of freedom from the 
value of \ 2 for the best-fit model. This estimate was compared with 
a covariance matrix built from multiple realizations of £ (it, n) us- 
ing Gaussian fields; consistent errors were obtained. We applied 
the simple method of Peacock et al. (2001) to the current data, 
keeping the assumed APM f (r), and obtained the marginalized re- 
sult /3 — 0.55 ± 0.075. The comparison with our best estimate of 
(3 = 0.49 ±0.09 indicates that the systematic errors in the previous 
analysis (from e.g. the assumed £[r]) were not important, but that 
the previous error bars were optimistic by about 20%. 

7.4 Peculiar velocities as a function of scale 

There has been much discussion in the literature on whether or not 
the pairwise peculiar velocity dispersion, a, is a function of pro- 
jected separation, a. Many authors have used A-body simulations 
to make predictions for what might be observed. Davis et al. (1985) 
found that the pairwise velocity dispersion of cold dark matter re- 
mains approximately constant on small scales, decreases by about 
20-30% on intermediate scales and is approximately constant again 
on large scales. Cen, Bahcall & Gramann (1994) found a similar 
overall behaviour as did Jenkins et al. (1998) whose results are plot- 
ted in the left panel of Fie. l23l as the solid line for a ACDM cosmol- 
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Figure 23. The variation of a with projected separation, a. Left panel: 2dFGRS data compared to some analytical models, as indicated in the legend. Right 
panel: 2dFGRS data compared to other redshift surveys and simulated catalogues, as indicated in the legend. The 2dFGRS data and simulated catalogue results 
use /3 to calculate the infall velocities whereas the other results assume a functional form (see discussion in Section lT4l . 



ogy. The dashed line is from Peacock & Smith (2000) who used the 
halo model to predict the peculiar velocities for the galaxy distribu- 
tion. Kauffmann et al. (1999) and Benson et al. (2000) used the GIF 
simulations combined with semi-analytic models of galaxy forma- 
tion and the galaxy predictions of Benson et al. (2000) are shown 
by the dotted line. These predictions generally assume erg = 0.9, 
but there is evidence that erg could be 10% lower than this (Spergel 
et al. 2003) and so the pairwise velocity dispersions implied would 
also be lower. 

Observationally, Jing et al. (1998) measured the pairwise ve- 
locity dispersion in the Las Campanas Redshift Survey and found 
no significant variation with scale. We note again that the errors 
for the LCRS ignore the effects of cosmic variance and are likely 
to be underestimates. Zehavi et al. (2002) used the SDSS data and 
found that a decreased with scale for a > 5 /i _1 Mpc. These ob- 
servations are plotted in the right panel of Fig. 1231 All these obser- 
vations have assumed a functional form for the infall velocities (or 
'streaming') and not used (3 directly. We have already shown that 
proper consideration of the infall parameter is vital in such stud- 
ies. Indeed, Zehavi et al. (2002) say that their estimates of a for 
a > 3 h~ 1 Mpc depend significantly on their choice of stream- 
ing model. This factor, along with a dependence of a on luminosity 
and galaxy type may help to explain the differences between the 
2dFGRS and SDSS results. 

The difference in results from Section lo"2l which measured the 
value of a at small a (570 km s~ ), and from using the £(cr, 7r) 
grid (506 km s _1 ), which measures an average value, hints that 
there may be such an dependence of a on a in the 2dFGRS data. 
We test for variations in a by repeating the fits described in Sec- 
tion l7.1l using a global /3, ro and ~/ r but allowing a to vary in each 
a slice. The results are shown in Fig. 1231 compared with the re- 
sults from other surveys, and numerical simulations as discussed 
above. The value of 506 km s _1 obtained from the 2-d fit for scales 



> 8 hT Mpc is close to the value at 8 h~ Mpc where most of the 
signal is coming from. The value of 570 km s _1 obtained from the 
Fourier transform technique agrees well with the results found for 
a < 1 /i _1 Mpc. The values of /3, ro and "/ r are essentially un- 
changed when fitting in this way. We note again that the effects of 
the infall must be properly taken into account in these measure- 
ments. We also note that we used our linear, power-law model on 
all scales, but we have seen that this is a reasonable approximation 
on non-linear scales. 

We see that the overall shape of the 2dFGRS results are fairly 
consistent with, though slightly flatter than the semi-analytic pre- 
dictions, but the amplitude is certainly a little different, which could 
be due to the value of erg used in the models, as discussed above. 
We also plot the mean of the mock catalogue results (solid line), 
and the results of a simulated catalogue (dashed line) of Yang et al. 
(2003, with as = 0.75) and these match the real data well. 



8 CONSTRAINING O m 

We take the value of (3 measured from the multi-parameter best fit 
to£(er,7r), 

0(L s ,z s ) = 0.49 ±0.09, (21) 

which is measured at the effective luminosity, L 3 , and redshift, z s , 
of our survey sample. In Section l2~T1 we quoted these values, which 
are the applicable mean values when using the J3 weighting and 
redshift cuts employed, as L 3 ~ 1.4L* and z 3 ~ 0.15. We also 
note here that if we adopt an Sl m = 1 geometry we find that (3 = 
0.55, within the quoted lcr errors. 
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8.1 Redshift effects 

The redshift distortion parameter can be written as, 

_ /(n m ,n A ,z) 



(22) 



where / = dlnD/dlna, and D is the linear fluctuation growth 
factor and a is the expansion factor. A good approximation for /, 
at all z, in a flat Universe, was given by Lahav et al. (1991), 



(2 - fi m - fi™)/140 « n° 



(23) 



and so to constrain f2 m from these results we need an estimate of b. 
There have been two recent papers describing such measurements. 

Verde et al. (2002) measured b(L s ,z s ) from an analysis of 
the bispectrum of 2dFGRS galaxies. Their results depend strongly 
on the pairwise peculiar velocity dispersion, a, assumed in their 
analysis. They used the result of Peacock et al. (2001), who found 
a = 385 km s , somewhat lower than our new value of « 
500 km s . To derive f2 m using these results would not therefore 
be consistent and so a new bispectrum analysis is in preparation. 

Lahav et al. (2002) combined the estimate of the 2dF power 
spectrum, P(k) (Percival et al. 2001), with pre- WMAP results from 
the CMB to obtain an estimate of b, but this value is also dependent 
on f2 m . Their likelihood contours 1 are reproduced in Fig. 1241 as the 
dashed lines. They also introduced a 'constant galaxy clustering' 
model for the evolution of b with z. Following these equations we 
can evolve our measured j3 to the present day and estimate 



(3{L s ,z = 0) = 0.45 ±0.08 



(24) 



and these contours are shown by the solid lines in Fig. [24] These 
are in good agreement with, and orthogonal to, those of Lahav et 
al. (2002). 

8.2 Luminosity effects 

We note that the above analysis is independent of luminosity as we 
examine everything at the effective luminosity of the survey, L s . 
From the correlation functions in different volume-limited samples 
of 2dFGRS galaxies, Norberg et al. (2001) found a luminosity de- 
pendence of clustering of the form (cf. Eqn. llO> . 



b/b* =0.85 + 0.15(L/i*) 



(25) 



which gives an estimate for the bias of the survey galaxies, b 3 = 
1.066* (using L = 1.4L*), where b* is the bias of L* galaxies. If 
this bias relation holds on the scales considered in this paper then 
P will be increased by the same factor of 1.06, 



0.52 ± 0.09 



(26) 



and evolving j3 in a 'constant galaxy clustering' model (Lahav et 
al. 2002) then, 



0(L*,z = 0) = 0.47 ±0j 



(27) 



which we choose as a fiducial point to allow comparisons with 
other surveys with different effective luminosities and redshifts. La- 
hav et al. (2002) obtained, f3(L*,z = 0) = 0.50 ± 0.06, in their 
combined 2dFGRS and CMB analysis, completely consistent with 
the our result. 



The bias parameter measured by Lahav et al. (2002) depends on r, the 
optical depth of reionisation, as b oc exp(— r). The plotted results do not 
include this effect, which could be significant, and this is discussed further 
in Section lol 
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Figure 24. Constraints on Sl m and b - Solid lines: Best-fit and Icr error 
contours on from this work, evolved to the present day (see Section l8lTl . 
Dashed lines: lcr and 2cr error contours from Lahav et al. (2002). Dotted 
lines: lcr constraints from WMAP (Spergel et al. 2003). 



8.3 Comparisons 

Percival et al. (2002) combined the 2dFGRS power spectrum with 
the pre- WMAP CMB data, assuming a flat cosmology and found 
Clm(z = 0) = 0.31 ± 0.06. These measurements of fi m are also 
consistent with a different estimation from the 2dFGRS and CMB 
(Efstathiou et al. 2002) and from combining the 2dFGRS with cos- 
mic shear measurements (Brown et al. 2003). 

Also plotted in Figure [24] is the recent result from the analy- 
sis of the WMAP satellite data. Spergel et al. (2003) found fi m = 
0.29 ± 0.07 using WMAP data alone, although there are degenera- 
cies with other parameters. It is clear that this is completely consis- 
tent with the other plotted contours. Spergel et al. (2003) also found 
that the epoch of reionisation, r = 0.17, which would reduce the 
value of b found by Lahav et al. (2002) by about 16%, still in good 
agreement with the results in this paper. 



9 SUMMARY 

In this paper we have measured the correlation function, and vari- 
ous related quantities using 2dFGRS galaxies. Our main results are 
summarised as follows: 

(i) The spherical average of £(rr, n) gives the redshift-space cor- 
relation function, £(s), from which we measure the redshift space 
clustering length, so = 6.82 ± 0.28 h~ Mpc. At large and small 
scales £(s) drops below a power law as expected, for instance, in 
the ACDM model. 

(ii) The projection of £(<r, n) along the n axis gives an estimate 
of the real-space correlation function, £(r), which on scales 0.1 < 
r < 12 /i _1 Mpc can be fit by a power-law (r/ro)~ lr with ro = 
5.05 ± 0.26 ft _1 Mpc, 7 r = 1.67 ± 0.03. At large scales, £(r) 
drops below a power-law as expected, for instance, in the ACDM 
model. 
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(iii) The ratio of real and redshift-space correlation functions on 
scales of 8 — 30 ft -1 Mpc reflects systematic infall velocities and 
leads to an estimate of j3 = 0.45 ± 0.14. The quadrupole moment 
of £(<r, 7r) on large scales gives f3 — 0.47^o 'i6- 

(iv) Comparing the projections of £(cr, n) along the n and a 
axes gives an estimate of the distribution of random pairwise pe- 
culiar velocities, f(v). We find that large-scale infall velocities af- 
fect the measurement of the distribution significantly and cannot be 
neglected. Using f3 — 0.49, we find that f(v) is well fit by an expo- 
nential with pairwise velocity dispersion, a = 570 ± 25 km s~ , 
at small a. 

(v) A multi-parameter fit to £(<j, 7r) simultaneously constrains 
the shape and amplitude of £(r) and both the velocity distortion 
effects parameterized by j3 and a. We find (3 = 0.49 ± 0.09 and 
a — 506 ± 52 km s - , using the Hubble Volume £(r) as input to 
the model. These results apply to galaxies with effective luminosity, 
L ~ 1.4L* and at an effective redshift, z s ~ 0.15. We also find 
that the best fit values of /3 and a are strongly correlated. 

(vi) We evolve our value for the infall parameter to the present 
day and critical luminosity and find (3{L = L* , z — 0) = 0.47 ± 
0.08. Our derived constraints on £7 m and b are consistent with a 
range of other recent analyses. 

Our results show that the clustering of 2dFGRS galaxies as a 
whole is well matched by a low density ACDM simulation with 
a non-linear local bias scheme based on the smoothed dark-matter 
density field. Nevertheless, there are features of the galaxy distri- 
bution which require more sophisticated models, for example the 
distribution of pairwise velocities and the dependence of galaxy 
clustering on luminosity or spectral type. The methods presented 
have also been used on sub-samples of the 2dFGRS, split by their 
spectral type (Madgwick et al. 2003). 
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APPENDIX A: COHERENT INFALL EQUATIONS 

Kaiser (1987) pointed out that the coherent infall velocities take a 
simple form in Fourier space, 

P s (k) = (l+f3vl) 2 Pr(k). (Al) 

Hamilton (1992) completed the translation of these results into real 
space, 

£V, tt) = [1 + Pid/dzf^ 2 )- 1 ] 2 ^), (A2) 
which reduces to 

S'(a,w) = S (s)P ( f i)+b(s)P 2 (ix) + U(s)P 4 (n) : (A3) 
where in general, 

Us) = ( X + f + y) &)> (A4) 

6W= f¥ + ^n KW-IWL (A5) 



3 7 



*(r) + ~*(r) - |€(r) 



(A6) 



and 

CM = 4 r i(r')r' 2 dr\ (A7) 
r Jo 

!('') = 4 r ^y 4 dr'. (A8) 
r Jo 

In the case of a power law form for f (r) these equations reduce to 
the form shown in Eans, IT2l to 1 141 In the case of non-power law 
forms for the real-space correlation function these integrals must 
be performed numerically. 

This paper has been typeset from a TpX/ KTpX file prepared by the 
author. 
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